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Transient Radio-frequency Ground Waves Over the 
Surface of a Finitely Conducting Plane Earth 

J. R. Johler 

The complete transient signal is reconstructed after propagation via the ground- wave 
mode over a finitely conducting plane earth in which the displacement currents are neglected. 
The theory is illustrated by computations that have been made on formulas derived by the 
method of the inverse Laplace transformation. The results of this analysis indicate that 
current sources with sinusoidal form in the time domain could be used to simulate spherics. 
The methods employed in this analysis can be also used to reconstruct propagated signals 
of pulsed radio-navigation systems. 



1. Introduction 

In recent years considerable interest has developed 
in various types of electromagnetic ground-wave 
transients that propagate over the surface of the 
earth. This interest has in large measure been 
stimulated by the use of pulse techniques for pre- 
cision radio-navigation systems and the ever presence 
of spherics that are propagated around the earth 
from various thunderstorms originating principally 
in the vicinity of the tropical land masses of the 
earth. In a previous paper by the author [1] x and 
in papers by Wait [2, 3], a theory of propagation of 
the ground-wave transient was developed as an 
inverse Laplace transformation. The basic mathe- 
matical problem is the evaluation of the Fourier 
integral, which describes the transient signal, E(t' , d), 
in space and time as follows: 

E{t'4)=~ r°° exp (iut)f s (a>)E(u,d)da>, (1) 

^ J -oo 

where / s (<o) is the source current transform, and 
E(a>,d) is the transfer characteristic of the propaga- 
tion medium. The theory was illustrated [1] by an 
evaluation of the sinusoidal part of the signal. Such 
examples describe the complete signal at the larger 
values of time, at low frequencies, and over highly 
conducting earth. It was shown that a signal 
generated by a sinusoidal current source, which 
commences abruptly at some time, £=0, and which 
is allowed to continue oscillating indefinitely, reaches 
a steady state condition, the amplitude and phase of 
which correspond precisely to the continuous-wave 
type of signal, i. e., the signal not interrupted in time. 
It is the purpose of this paper to extend the theory and 
reconstruct the complete signal by the addition of 
certain nonsinusoidal terms, which appear as a result 
of the inverse Laplace transformation and which 
describe the early time behavior of the signal. 

2. Theory 

The transient signal, E(t' , d), may be represented 
over a finitely conducting plane earth in which the 



1 Figures in brackets indicate the literature references at the end of this paper. 



displacement currents are negligible, for a vertically 
polarized source of the form, exp(— erf) cos 
co c 2(0<£< oo ), at the local time, l' } as follows [1]: 



E(t,'d) = vC exp {—vt') < - 1— vJOci exp0 2 a) 

[erfc(-^)-erfc(^-^)]+^~^ 2 } 

+^^i^y°[M (2) 



where 
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and where 

d=t\\e distance along the surface of the earth, 

in meters, 
C=2(10~ 7 )/d, where IqI=1 ampere-meter [4], 
r/i=the index of refraction of air at the surface. 

ih ~ 1.000338, 
o-=the conductivity of the ground, in mhos/ 

meter, 
/loathe permeability of space, £i =47r (10~ 7 ) 

henry /meter, 
c=the speed of light, c-2.997925 (10 8 ) meters/ 

second, 
a=(rjid)/c; t—a=t', the local time, and 
p 1 r=the numerical distance, p^aco 2 , neglecting 

displacement currents in the earth. 

Also, the transient signal, E"(t' y d), for a source, exp 
(—Cif) cosco c £(0<£<r 2 ), i- e -> a source interrupted at a 
time, T 2 , may be represented as follows: 

E" (t',d) = E(JL' 9 d)-exp {-vT 2 ) E(t",d), (4) 
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t"=t'-T 2 



(t'>T 2 ). 



(5) 



The continuous wave or the wave not interrupted in 
time is the special case: £'>>0, Ci=0, T 2 = °° . The 
early values of time, £'>0, in this case describe the 
development of the steady state. 

The source function, F s (t), may be brought into 
existence less abruptly as follows: 



JF(0=exp(-^-exp(-#) 
Z=c 2 +iw C ' 



(6) 
(7) 



v has been previously described (2a), and c 2 is assigned 
a large positive value (e 2 >>e!). The transform, 
/(co), (1) may be written 

/(«) =/,(*,) £(«), (8) 

or more generally 

f(u)=E(v)V Stl (a>)+f 8t2 (<*)+ . . .], (9) 

where / s ,i(oj), / s . 2 (w), . . . are source functions that 
have been superposed to determine some complicated 
waveform. The solution of the space-time function, 
E'(t',d), for this source (6) then merely requires the 
sum of two waves calculated as previously de- 
scribed (2), 



E'=E'(fi',d)=E,(t',d)-Et(t',d). 



(10) 



Figures 4, 5, and 6 illustrate this signal. 

It is interesting to note the space-time function, 
E(t',d), for the special case, co c =0, or the case of a 
step function source, u(t), exponentially damped, exp 
(-erf) u(t)* 

E(t',d)=CiCexv (— crf'K — 1— CiVraexp (e?a) 
[erfc (-^)-erfc(^- Cl V^)]+^-^} 
+[ 2 4+c 1 ]^exp[=^]+c 1 c[^} (11) 



In this case, the imaginary part of the function, 
E(t' y d), vanished. 

The source has been assigned the form exp (—erf) 
cos wrf or without damping (ci=0), cos corf. The sine 
source, exp (—erf) sin corf merely involves the forma- 
tion of the product — i E(t',d) or the real part be- 
comes the positive imaginary +7m E(t',d). The sine 
source at zero frequency (11) therefore vanishes. It 
can be further concluded that the step function 
response of the ground-wave propagation medium 
for a step function current source, u(t), is the particu- 
lar case, C!=0, or v=0 (2), which is simply, 



E(t' ? d)=fcexp 
Za 



L4a J 



(12) 



3. Computation 

The detailed reconstruction of the complete tran- 
sient signal over a finitely conducting plane earth, 
neglecting displacement currents in the earth, is 
shown in figs. 1 to 6. The complete signal, E= 
E(t', d) (2), is presented as three terms, 

E=E(t',d) =E Q (f',d) +Ei (?,d) +E 2 (t',d) . (13) 

The first term, E =E (t', d), is called the sinusoid 
because it is multiplied by the time harmonic func- 
tion exp(— vt') (2). The remaining terms are 

El = E^4) = vC^ew[^y^y (14) 

(15) 



V 



=t Cex ^[-£ r __ 



E'=E'(t',d) K (figs. 4, 5, 6) describe a signal the 
source of which has been brought into existance less 
abruptly (6). 

Since the complete signal is very complicated, this 
presentation (13) illustrates the influence of the 
various parts of the signal as space-time and fre- 
quency parameters are varied. Thus, for example, 
the term E 2 degenerates into a "spike" close to the 
origin of time, t', as the frequency is decreased (fig. 
6), and vanishes for the less abrupt source (6) such 
that the terms Eq+E x suffice to describe the com- 
plete signal. This is consistent with the limiting 
condition of very small numerical distance (<r= oo) or 
an attenuation factor of unity. Thus, the trans- 
form of the Norton surface wave at considerable 
distance from the source may be written 



E(s) = l ( v =0 y a=0) y 

and the space-time function, E{t',d), is 
E(t',d)=^- 1 E(s)=8(t'), 



(16) 



(17) 



2 The step function, u(t), is defined as follows: u(t) = l for t>0; u(t)=0fovt<0. 



where 8(f) is the unit impulse function at the origin 
of time. Therefore, the transition to infinite con- 
ductivity illustrates the development of the Dirac 
impulse function, 8(t'). 

The reconstruction of a radio-navigation type of 
pulse is illustrated (fig. 7) with the "classical" sine- 
squared current source. The transient is allowed to 
reach the steady state because the "leading edge" 
of the pulse is of primary interest. It is common 
practice to interrupt radio-navigation pulses at a 
time, T 2y (4, 5) or apply exponential damping after 
some time, T 2 , has elapsed. The source is rede- 
fined as follows: 

Ee F s (t)=sm 2 a> p t sin corf (0<rf< oo ), (18) 



i % 

F,{t)=2 exp (— vt) — | exp (— Vl t) 



-exp (—v 2 t), 
(19) 
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Figure 1. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 
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Figure 2. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 
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Figure 3. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 
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Figure 4. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 



283 



— 2 



<u o 



Complete Signal, E( Amplitude) 


'' 


VERTICAL ELECTRIC DIPOLE SOURCE 
I i - 1 Ampere -meter 
lu f e = 100 kc 
, ,N k d = 10 Miles 
f \ or * 0.005 Mho per Meter 
1 \ p. • 0.0188 
| V c, =2.5 (10 5 ) 
^ , ^^-^Plane Earth Theory 
En + E f {Amnlltiirtp)~~~ — 


Note: E = E + E,+E 2 
E'=E + Ef 

c 2 =a5(io s i 

Where 
Re F s (t)= [expt-ctl-expl-CztJlcosa-ct 


1 I ....... , ,.,..,, , 



-Complete Signal, E( Phase) 



KIQ-5) 


1 "'_ 




Y~ Complete Signal, E 


DUO" 6 , 


] 




jl / \ ^ _ 




^^i- Complete Signal, E' 




I I | ''■■■'■■■ ■ 1 I 1 1 I L.1 ! 1 1 1 1 1 1 , 1 I.I 



8 I 




Time, t', Microseconds 

Figure 5. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 
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Figure 6. Amplitude, phase, instantaneous signal, and cor- 
responding source current for propagated surface-wave 
transient. 
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Figure 7. Amplitude, phase, instantaneous signal, and cor- 
responding sine squared source current for propagated surface- 
wave transient. 



where 



v=Ci+io: c (as before), 

vi=c l +i(u c +2u P ), 
v 2 =Ci+i(a> c —2w p ). 



(19a) 
(19b) 
(19c) 



The exponential damping factor, c u may be assigned 
a finite value or zero. The space-time function, 
E(t',d), may therefore be written as the sum of three 
oscillatory waves; 

E(t\d)=^E v (t'4)-lE yl (t\d)^E v2 (t\d). (20) 

Figure 7 describes the propagated wave (20). Each 
term is calculated as described previously (2) 
(ci=0). The dispersion of the pulse at early times 
is a result of the transient behavior of the ground 
wave. The dispersion of the pulse at later times is a 
result of the operation of the propagation medium on 
the widely separated frequency spectra that develop 
after considerable time, ¥ , has elapsed. It is in- 
teresting to suggest that the amplitude and phase 
could simulate the output of amplitude and phase 
envelope detectors of a radio receiver. However, 
the finite bandwidth of the receiver, / r (co), should 
be introduced into the transform, /(w), (9) as follows: 

/(«)=£(co)[/, il («)+/, i2 (co)+ . . . ] 



F/r,l(")/r,2(w) 



L (21) 



where /r,j(co), ,f r)2 (co), . . . are the transfer char- 
acteristics of the frequency selective networks in 
the receiver circuit. 
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4. Conclusion 

The propagated transient signal over a finitely 
conducting plane earth can be reconstructed in de- 
tail, for current sources of complicated form, for all 
values of local time, by application of the inverse 
Laplace transformation. The most important merit 
of this analysis is the insight gained into the details 
of the propagation mechanism. It should be 
emphasized that the reconstruction of the complete 
signals at times less than approximately 0.1 micro- 
second requires the introduction of displacement cur- 
rents and earth curvature, which have been neglected 
in this paper [3]. The method does however serve 
as a check on the tedious numerical methods of inte- 
gration of the Fourier integral, the direct evaluation 
of which seems to be the obvious method of intro- 
ducing into the transient solution of the ground wave 
the effect of the earth curvature and the effect of the 
displacement currents in the earth. 



The mathematical formulas were programed for 
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Cambridge, Mass., and the Radio Navigation project 
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